Nonequilibrium quantum dynamics of a charge carrier doped into a Mott insulator 
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We study real-time dynamics of a charge carrier introduced into undoped Mott insulator propa- 
gating under a constant electric field F on the t-J ladder and square lattice. We calculate quasis- 
tationary current. In both systems adiabatic regime is observed followed by the positive differential 
resistivity (PDR) at moderate fields where carrier mobility is determined. Quantitative differences 
between ladder and 2-dimensional (2D) system emerge when at large fields both systems enter nega- 
tive differential resistivity (NDR) regime. In the ladder system Bloch-like oscillations prevail, while 
in 2D the current remains finite, proportional to l/F. The crossover between PDR and NDR regime 
in 2D is accompanied by a change of the spatial structure of the propagating spin polaron. 

PACS numbers: 71.27, 71.27.+a, 72. 10.Bg 
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Introduction. — The real-time response of interacting 
many-body quantum systems remains in many aspects 
an unexplored field that has recently attracted signifi- 
cant attention. Understanding of the time-dependent 
phenomena becomes vital for various branches of physics 
including condensed matter nanostructures [2| and 
optical lattice systems 0,0] ■ Since only very few cases are 
exactly solvable, the vast majority of unbiased results has 
been obtained from numerical approaches like exact diag- 
onalization (ED) 0, @, time-dependent density matrix 
renormalization group Q or nonequilibrium dynamical 
mean-field theory Among others, the electric-field 
induced breakdown of the Mott insulator (MI) la Ifil- 
ll| , nonequilibrium transport in nanostructures [121 1 1 31 ] 
and the relaxation of correlated systems after the pho- 
toexcitations [13, 15|. represent the well-studied exam- 



ples of nonequilibrium phenomena which are important 
both for fundamental understanding of strongly corre- 
lated systems as well as for their potential applications. 

Most of theoretical studies so far considered the break- 
down of undoped MI, when threshold value of the electric 
field exceeds experimental value |9|1 by a few orders of 
magnitude (see discussion in Ref. It indicates that 

other transport mechanism becomes active at energies 
much lower than the Mott-Hubbard gap. In this Letter 
we inveqstigate a nonequilibrium response of a charge 
carrier doped into the insulator and driven by an uni- 
form electric field F. Understanding of this subject on 
one hand widens our knowledge of a charge carrier doped 
into the antiferromagnetic ( AFM) background H3] , 
on the other, it represents a fundamental problem of a 
quantum particle moving in a dissipative medium [lil ] . 

Having in mind strongly correlated systems, we inves- 
tigate the t-J model where the particle driven by a con- 
stant electric field dissipates the energy by inelastic scat- 
tering on spin degrees of freedom. We use numerical ap- 
proaches to treat t-J model at zero temperature on two 
different system geometries, i.e., a ladder with periodic 
boundary conditions and an infinite 2-dimensional (2D) 



square lattice. 

The most important finding of this Letter concerns the 
current-field characteristic, where the linear part with a 
well defined mobility is followed by a strongly nonlin- 
ear one with a negative differential resistivity. Numeri- 
cal results reveal also adiabatic evolution for very weak 
fields (below the linear regime) and strong Bloch oscilla- 
tions (BO) in the opposite limit of large F. These qual- 
itative results hold for both the geometries. The most 
prominent difference between the ladder and the 2D case 
emerges at large field where in the latter system due to 
different topology, allowing for transverse charge carrier 
motion, BO remain damped and the steady current de- 
creases with field as l/F. 

Model — We consider a charge carrier within the t-J 
model threaded by a time dependent magnetic flux: 
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where Cj i(T = Cj i(T (l — n- h — a ) is a projected fermion op- 
erator and (lj) denote nearest neighbors. The constant 
electric field F is switched on at time t = and is mea- 
sured in units of [to/eoa], where eo is the unit charge and 
a is the lattice distance. We set to = h = bq = a = 1. 
For the ladder the charge current operator for t > reads 
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where F acts along the ladder's leg and <^>ij(i) equals —Ft 
and for hopping in x- and y-direction, respectively. 

At this stage it is instructive to recall a simple relation 
between the current j(t) = (I(t)) and the total energy 
E(t) = (H(t)) I9j], E(t) = FI(t), which allows one to 
calculate the steady component of j(t) 
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where AE(t) = E(t) - E(0). 
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Ladder. — The real-time response of a ladder with L 
rungs is studied in the full Hilbcrt space by means of ED. 
Applying the Lanczos technique we have determined the 
initial ground state = 0)) at p = (tt/2, 0) |2(J. The 
time evolution of the initial state is calculated by step- 
vise change of the flux 4»ij{t) in small time increments 
5t <S^l 7 employing at each step Lanczos basis generating 
the evolution \V(t-6t)) -> \$(t)) [3[21|. 

Fig. [TJa) demonstrates how the total energy changes 
in time. One can identify three regimes/limits of F: the 
adiabatic regime (AR) for F — > 0, the Bloch-oscillations 
regime (BR) for F>1, and the dissipative regime (DR) 
for intermediate F. In the adiabatic limit, E(t) follows 
the ground state dispersion -Eo(p) with p = Po + Ftx. 
In the ladder E'o(p) is separated from excited states by a 
finite gap, related to spin gap. The dispersion is a quasi- 
particle (QP) one, i.e. it follows the behavior of a single 
carrier with the periodicity Ap x = lix. As a result E{t) 
oscillates with the Bloch frequency ojb = 2ir/tB = F also 
in AR. Due to finite gap, qualitatively similar behavior 
remains even for \F\ > 0. We expect that the crossover 
from AR to DR is connected with a transition through 
the gap resembling the Landau-Zener transition through 
the Mott-Hubbard gap @. Hence, the threshold F 
should be determined by J. 

Fig. HJa) shows that E(t) is monotonic only in DR 
where j{t) doesn't change sign. Other regimes (AR, BR) 
are characterized by strong current oscillations. Fig. 
[Ha) together with Eq. ([3]) allow one to estimate j 
as a function of the electric field. Clearly j is maxi- 
mal in DR, therefore there exists a corresponding "op- 
timal" F. Substantial differences between the regimes 
show up also in the kinetic energy related with the move- 
ment along the i-direction, Ef. x — (Hkx) where Hk x = 
_ *o I2j.<j (c~ lFt c] +i)(r Cj> + H.c.J . We have found that 
Ekx(t) oscillates in AR and BR. This quantity is always 
negative in AR, whereas in BR it takes on negative as 
well as positive values. Therefore, BR resembles the BO 
of nonintcracting particles in that both j(t) and Ek x (t) 
change sign. 

Nonzero j in DR implies a steady growth of energy 
due to the Joule heating (although we are not dealing 
with well defined thermalization). This effect is limited 
by the system size since a carrier driven by constant F 
propagates many times around the ladder with finite L, 
thus steadily excites the spin background and increases 
its effective temperature. Fig. QJa) shows that AE(t) 
deviates considerably between different L at time t\ ~ 
21b- As there is a single charge carrier, the quantity 
r(t) = f* drj(r) = AE(t)/F [see Eq. ©] can be viewed 
as a distance traveled by a carrier within the time-interval 
(0, t). For noninteracting case AE(t) < 4 and r(t) < 
i/F, which is the Stark localization length [22j]. Here, 
we have found that r(ti) ~ L confirming that the finite 
size effects originate from heating of the spin background 
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FIG. 1: (Color online) Ladder with L rungs: (a) increase of 
energy AE(t) in different regimes; (b) real-time current j(t) 
normalized by kinetic energy a(t) = j(t) / (F\Ekx(t)\) (note 
the steady ^-independent ratio); (c) kinetic energy E^ x vs. 
tits, (inset) and t/Lts (main); (d) steady current j vs. F. 
Dotted line in (a) shows AE(t) estimated from the mobility 
fi shown in (d) as linear fits j — fiF. J = 0.6 is used in (a)-(c). 



when carrier repeatedly encircles the ladder. 

However, even for t > ti one can follow the behav- 
ior in a controlled way. Previous analysis of ID system 
of spinless fcrmions has shown that heating can be ac- 
counted for by a renormalization of the kinetic energy, as 
a sum rule for the optical conductivity <j(u>) Re- 
sults in Fig. [ljb) show that the same reasoning can 
be applied to the t—J ladder. In particular, the ratio 
a (t) = j(t)/(F\Ekx(t)\) reaches already at short t/ts 
nearly a constant value that is independent of F and 
L. Therefore, the original problem concerning the evalu- 
ation of j can be reduced to finding E£ x = Ek x (t ~* 00 ) 
in the limit L — > oo. 

A straightforward expectation is that for n-times 
larger system [L — > nL) it takes n times longer to reach 
the same average temperature of the spin background. It 
explains the E^ x oc Ek x (t/L) scaling visible for t > t\ in 
Fig. Hljc). The initial value of E^ x in this time-domain 
represents the upper bound on E£ x when L — > oo. On the 
other hand, for t < t\ results for Ek x (t) merge without 
any rescaling of time [see the inset in Fig. QJc)]. Hence 
the final value of E^ x in the latter time-domain poses 
the lower bound on E£ x . These bounds on E£ x together 
with the ratio a(t) allow one to extract j (see Fig. Old)). 
Within presented range of F, in the DR nonlinear effects 
are weak and one can easily estimate the hole mobility: 
H ~ 1.7 for J = 0.3 and /i ~ 1.3 for J = 0.6. However, 
significant nonlinear effects have to show up for larger 
fields since j is maximal for a finite F. This observation 
together with the value of are the main results for the 
ladder system. 
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2D square lattice. — Since in the ground state a single 
carrier in 2D lattice carries momentum ko = (7r/2,7r/2), 
we set </>ij(i) = —Ft/y/2 for x- as well as y-direction, 
and ujb = F/y/2. Accordingly, we also define and cal- 
culate the modified current j(t), Eq. @, along the di- 
agonal. We employ an ED method at defined over a 
limited functional space (EDLFS) which describes prop- 
erties of a carrier doped into planar ordered AFM [171 ]. 
One starts from a translationally invariant state of a car- 
rier in the Neel background \<f>o) = Ck |Neel). The kinetic 
part Hk as well as the off-diagonal spin-flip part Hj of 
the Hamiltonian (Q} are applied up to 2V/, times generat- 
ing the basis vectors: {\4>\ Lh )} = [H k (F = 0) + Hj] nh \<p ) 
for rih = 0, ...,2V/,. The ground state \^(t = 0)) and its 
time-evolution under electric field are calculated within 
the limited functional space in the same way as previously 
for the ladder within the full Hilbert space. The advan- 
tage of EDLFS over the standard ED follows from sys- 
tematic generation of selected states which contain spin 
excitations in the vicinity of the carrier. It enables in- 
vestigation of the dynamics of large systems, which are 
far beyond the reach of ED. Since 2V/, determines the ac- 
cessible energies (spin excitations), this quantity poses 
limits on the maximal propagation time characterized by 
a steady growth of energy. We set 2V/, = 14 while smaller 
values of 2V/, are used to establish the time-window where 
results are independent of 2V/, (see inset in Fig [3k) . 

AR at small F = 0.1 is clearly seen from Figs. [Ha) 
and (b) that show AE(t) and j(t) after the field has 
been switched on at t = 0. Both quantities are consis- 
tent with the adiabatic propagation along the QP band 
with a period t p = £b/2, consistent with the AFM long 
range order. In the DR the carrier starts to move due 
to the constant field, it emits magnons and consequently, 
after propagation through the transient regime t/ts <C 1, 
it develops a finite average velocity as seen for F = 0.6 
and 1.4. DR is characterized by a constant current and 
a linear increase of the total energy of the system. To 
calculate the current we make use of Eq. where the 
linear increase of AE(t) provides the value of j. Indeed, 
linear dotted-dashed fits in Fig. [2]^a) indicate that the 
system has already reached the quasi stationary state. 
We plot the extracted values of j as dotted-dashed hor- 
izontal lines in Fig. &b). Similar approach has recently 
been applied to the problem of a driven Holstcin polaron 
[lij . At large field, F > 3, jit) is consistent with a 
damped BO with t p < is, signaling the onset of BR. 

Fig. [2(c) displays j—F characteristics in DR for J = 
0.3 and J = 0.6. Each point has been calculated from 
Eq. ([3]) and then compared with j(t), as demonstrated 
in Figs. [2(a) and (b). The most remarkable property 
is a peak at Fq dividing DR into two sub- regimes, i.e., 
the regime of positive differential resistivity (PDR) for 
F < Fo and the regime of negative differential resistiv- 
ity (NDR) for F > Fq. The value of the crossover field 
Fq scales with the exchange energy J. From curves in 
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FIG. 2: (Color online) 2D lattice: (a) increase of energy 
AE(t) for J = 0.3, various F (main) and N h = 12, 13, 14 (in- 
set); (b) real-time current j(t), j from Eq. [3] (dotted-dashed 
lines); (c) steady current j vs. F fitted by j = /j,F (dashed 
lines) and j = b/F (solid lines) in PDR and NDR regimes, 
respectively. 



Fig. [2(c) we find Fo ~ 2.3J. This is consistent with an 
intuitive expectation that the onset of the NDR regime 
emerges when the energy gained by a single hop exceeds 
the maximal energy of one-magnon excitation. Conse- 
quently, the real-space propagation of a carrier exhibits 
qualitatively different behavior in both regimes. This 
is illustrated in Fig. [3] which displays spin deformation 
function C(r) around the carrier at different times. We 
define C(r) = | S^ el - S? +I |) with 5/ ' N ' cl = ±±. 

In the PDR regime and for F = 0.6 spin excitations 
predominantly emerge behind the traveling carrier indi- 
cating that the average charge velocity v c is larger than 
magnon velocity v s = v2J, i.e. v c = j > v s (see the 
upper panel of Fig. [3]). A more complex pattern charac- 
terizes the NDR regime where at F = 3, v c < v s (lower 
panel of Fig. [3|) and enhanced spin excitations develop 
also in front of the carrier. Moreover, C(r) displays a 
distinct transverse orientation with respect to the field 
direction. This signals an increased transverse and longi- 
tudinal carrier oscillation that serves to release the excess 
energy to the unperturbed spin background. 

In the PDR regime (dashed lines in Fig. [2(c)) linear 
fits of the current increase provide an estimate for the 
carrier mobility, which yields fx ~ 1.3 for J = 0.3 and 
(i ~ 0.8 for J = 0.6. These values are in agreement with 
ED calculations using the linear-response theory [23| and 
indicate that carrier mobility increases with growing cor- 
relations. In the NDR regime a scaling j ~ \/F is found, 
Fig. 12(c). This is consistent with incoherent hopping be- 
tween Stark states appropriate for dispersive boson exci- 
tations [24]], in contrast to coherent hopping in the case 
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FIG. 3: (Color online) 2D lattice with J = 0.3. Snapshots 
of the spin deformation C(r) in the vicinity of hole [placed at 
(0,0)], taken at different times t/ts- Filled arrows indicate 
direction of the electric field. 



of dispersionless phonons leading to j oc 1/ \/F [l8| • 

Discussion. — At small F the ladder and the 2D system 
exhibit adiabatic propagation with j = 0. The validity of 
AR is closely connected with the stability of the lowest 
QP band. While this is more evident for the ladder (due 
to a gap), the gapless (acoustic) magnons in 2D lattice 
require more care. Nevertheless, our results show that 
the condition for Cerenkov radiation of magnons v(k) > 
v s = \[2J is for J > 0.3 never fulfilled. In part, because 
the QP bandwidth is limited by J and not by t. For 
J = 0.3 we obtain maximal charge velocity i> max = 1.36 J 
while for J = 0.6, w m ax = 1.10J. Moreover, the QP 
wei ght remains finite throughout the whole Brillouin zone 
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With increasing F both systems enter DR where the 
quasistationary current is proportional to F, leading to 
well defined mobility being even quantitatively close for 
both cases. The mechanism of dissipation is clearly the 
emission of spin excitations. Still the dissipation is due 
to topological difference (in particular due to the pos- 
sibility of spin perturbation transverse to the direction 
of the carrier propagation) more efficient within the 2D 
lattice leading to stronger damping of the BO and finite 
j cx l/F. Consequently, quasistationary current on the 
square lattice occurs up to much larger F than in the 
ladder where in the same regime strong BO with j ~ 
prevail. 



Surprisingly, the undoped [1 1| and lightly doped MI 
show the same sequence of the field-regimes: j is zero 
or exponentially small for sufficiently weak F (AR in the 
present case); for larger F the linear j — F dependence 
is restored (DR); finally, for very large F the response 
is dominated by the BO (BR). However, in doped and 
undoped MI these regimes occur at exceedingly different 
fields. In particular, the DMFT studies of the strong-?/ 
limit 111 reveal that the threshold field for the dielectric 



of magnitude above the field discussed here. 
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